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Abstract 

The "extended Overhauser model" [Overhauser, A.W. Can. J. Phys. 
1995, 73, 683] for the calculation of the spherically and system-averaged 
pair density (APD) has been recently combined with the Kohn-Sham 
equations to yield realistic APD and correlation energies. In this work 
we test this approach in the high-density (weakly-correlated) limit of the 
He isoelectronic series and of the Hooke's atom isoelectronic series. Unlike 
many of the commonly used energy functionals, the Overhauser approach 
yields accurate correlation energies for both series. 

1 Introduction 

Kohn-Sham (KS) Density Functional Theory 1-3 (DFT) is nowadays one of the 
most popular methods for electronic structure calculations both in chemistry 
and solid-state physics, thanks to its combination of low computational cost 
and reasonable performances. The accuracy of a KS-DFT result is limited by 
the approximate nature of the exchange-correlation energy density functional 
E xc [n]. Simple approximations (local-density approximation and generalized 
gradient corrections) for E xc [n] provide practical estimates of thermodynam- 
ical, structural and spectroscopic properties of atoms, molecules and solids. 
However, with the current approximations, KS-DFT is still lacking in several 
aspects, in particular it fails to handle near-degeneracy correlation effects (rear- 
rangement of electrons within partially filled shells) and to recover long-range 
van der Waals interaction energies. The inaccuracy of KS-DFT stems from 
our lack of knowledge of E xc [n], and much effort is put nowadays in finding 
new approximations to this term (for recent reviews, see, e.g., Refs. 2-4). A 
popular trend in the development of new KS E xc [n] is the use of the exact ex- 
change functional E x [n] (in terms of the KS orbitals) , and thus the search for 
an approximate, compatible, correlation functional E c [n]. 
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An exact expression for E c [n] is the coupling-constant integral 5 ' 6 

E c [n} = J dXj^ dr l2 Anr\ 2 j^r l2 ) ^ , (1) 

where the interaction between the electrons is adiabatically turned on from 
u> A=0 (ri2) = to the Coulomb repulsion w A=Aphya (ri2) = l/ri2 by varying a 
real parameter A (typical examples are w x (r 12) — A/ri2, with A p h ys = 1, or 
w x (ri2) — erf(Ari2)/ri2, with A p h ys = 00). The one-electron density n(r) is 
(ideally) kept independent of A and equal to the one of the physical system by 
means of a suitable external potential v x (r). In Eq. ([T]) the correlation part of 
the spherically and system-averaged pair density (intracule density) f x (r 12) is 
defined as follows. For each A, take the square of the many-electron wavefunction 
*f? x ground-state of the hamiltonian H A , 

N _ 2 iV JV 

^ A = -E^+2 £ ^(Iri-rjD + XXfc), (2) 

i—1 i^j—l i—1 

and integrate it over all variables but the scalar electron-electron distance r 12 = 
ki -r 2 |, 

f\r l2 ) = N{N 2 l) £ J\*Hr u ,K,r 3 ,...,r N )\ 2 ^dKdr 3 ...dr N , (3) 

(Tl ...(TN 

whereR= (ri+r2)/2. The correlation part / A (ri 2 ) is then defined as / A (r 12) = 
f x (r 12) — /ks('*12), where the intracule density of the KS system is /ks(?"12) = 
/ A=0 (ri 2 ) (and yields the Hartree plus the exchange energy). 

The traditional DFT approach to the construction of approximate E c [n] 
is based on the idea of universality. For example, the familiar local-density 
approximation (LDA) consists in transfering, in each point of space, the pair 
density from the uniform electron gas to obtain an approximation for / A (r 12) in 
Eq. |T]). In a couple of recent papers, 7-9 we have started to explore a different 
way of constructing 2? c [ti], based on an "average pair density functional the- 
ory" (APDFT), which was inspired by the seminal work of Overhauser 10 and 
its subsequent extensions. 11-13 In this approach, we solve a set of radial (one- 
dimensional) Schrodinger-like equations that give, in principle, the exact f x {r\ 2 ) 
along the DFT adiabatic connection. In practice, this formalism contains an 
unknown effective electron-electron interaction that needs to be approximated. 
The APDFT equations must be solved for each system, and combined self- 
consistently with the KS equations. 9 Preliminary applications of this approach, 
combined with a simple approximation 7 for the effective electron-electron inter- 
action that enters in the formalism, gave accurate intracule densities f(r 12) and 
correlation energies E c [n] for the He isoelectronic series. 7,9 

Katriel et al. 14 have recently tested most of the currently available correla- 
tion energy functionals in the high-density (weakly-correlated) limit of the He 
and of the Hooke's atom isoelectronic series, finding that, while several function- 
als are accurate for the He sequence, none is satisfactory for the Hooke's atom 
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series. Motivated by their findings, in this work we compute the correlation 
energy and the intracule density in the high-density limit of the two series via 
the APDFT approach, 7-9 finding accurate results in both cases. 

The paper is organized as follows. In the next Sec. Owe recall the basic equa- 
tions that define the high-density limit of the He and Hooke's atom sequences, 
to which we apply, in Sec. [3] and 01 the formalism of Refs. 7-9 to compute the 
intracule density and the correlation energy. In Sec.[5]we also analyze the failure 
of LDA in this limit from the point of view of f(r 12)- The last Sec. |5]is devoted 
to conclusions. 



2 The high-density limit of the He and Hooke's 
atom isoelectronic series 

The two hamiltonians analyzed in this paper read 

H = -^-^+v(r 1 )+v(r 2 ) + —, (4) 
2 2 ri2 

/ s f — — He series 

V ( T ) = \ r 

\ \kr 2 Hooke's atom series. 

We are interested in the high-density (weakly-correlated) limit, which corre- 
sponds to Z — > 00 and k — * 00. By switching to scaled coordinates s = r/a, 
with a = Z^ 1 (He series) and a = fc -1 / 4 (Hooke's series), both hamiltonians 
have the form 

H = A (-^T -^T + ^ + C ^ + + offi), (5) 

cr \ 2 2 S12 / cr 

where v(s) = —1/s for the He series, and v(s) = s 2 /2 for the Hooke's atom 
series. We thus study pertubatively the system described by Hq + otH\. 

The order zero of the one-electron density n(r) and of the intracule density 
/(ri2), in scaled units, is simply 

mw , f ^e" 28 (He) 

ni )(S) = { ^ (Hooke) <«> 

fi o )(q , J ^(3 + 6 Sl2 + 4 Sl2 )e^- (He) 

1 (Sl2) = \ j^me-^ (Hooke) ^ 

These functions are correctly normalized, so that if we switch back to coordi- 
nates r we have n^(r) = a~ 3 n^°'(s — a -1 r), etc. 

The first-order correction to the scaled density, n(s) — (s) + a n^ 1 ' (s) + ... 
is given by 

n^(s) = 2n^(s) X (s), (8) 

where 15 

/% 23c 3 31 ■ c 5 3 __ . , „ . 3,,, 
X(s) =--- — - - 7 + -— — + -s + -E l( -2,s) - - ln( S ), (9) 
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for the He isoelectronic series, with 7 = 0.577216.., 



and 16 
X(s) 



Ei(-x) = 
erf(s) V2(l + ln2) 



1 



e^^erfc 



x 

72 



sV^ Jo 

/V2 



-dt, 



dx 



(10) 



( e -(*-) 2 _ e -(^) 2 ) 



*e t2 erfc(t) 



(11) 



for the Hooke's atom isoelectronic series. 

By definition, the Kohn-Sham hamiltonian describes a non-interacting sys- 
tem that has the same density of the physical, interacting, system. Thus, the 
first-order change in the electron density of Eq. ((5|) corresponds to a first-order 
change in the KS system. Therefore, we write the scaled intracule /(S12) up to 
orders a as 



/M = / (0) (*i 2 ) + a \f§£M + fP(s 12 )] + O (a 2 ) 



(12) 



where we have separated the first-order correction into a Kohn-Sham part and 
a correlation part. The KS part /4g is entirely determined by the first-order 
density of Eq. ©, 



dsn^(s + s 12 )n^\s)x(s), 



(13) 



and is reported in Appendix [A] in analytic form for the He isoelectronic series, 
while is obtained numerically for the case of the Hooke's series 

The total first-order intracule f^ 1 ' 
case of the Hooke's series, 16 



fits + /<= i s known analytically in the 



/ (1) (.si 2 ) = 



2e"^/4 
(2tt) 3 / 4 



1- 



1 + In 2 1 



2/2 



S12 



•S12 



— -l^-erfc(^)W2 



V2 



H2/V2 



e* erfc(t)dt 
(14) 



3 Effective equations for f(ru) in the high-density 
limit 

3.1 Formalism 

We are interested in calculating /j and the corresponding second-order corre- 
(2) 

lation energy Ec with the method of Refs. 7-9, in which the intracule density 
/(ri2) of the physical system is obtained from a set of effective equations, which 
for two-electron systems reduce to 

[-V 2 i2 + w eS {r 12 )} ij{r l2 ) = e V(ri 2 ), (15) 
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with /(r 12 ) = \ip{ri 2 )\ 2 . Equation can be derived by considering 8 ' 9 a set 
of Hamiltonians characterized by a real parameter £, 



iv 

E 



N 

E 



^(|r, - r,-!) +££i; ne (r i ), /«(r 12 ) = /(r 12 ) V£ 

i=l 

that describe a set of systems in which the external potential is turned off as 
£ — > 0, and the intracule density is kept fixed, equal to the one of the physical 
system, by means of a suitable electron-electron interaction w^(r 12 ). In the case 
N = 2, when £ = we have a translationally- invariant system (the center-of- 
mass degree of freedom is described by a plane wave) of two fermions in a relative 
bound state (similar to the case of positronium, but with a different interaction). 
This relative bound state is such that the square of the wavefunction for the 
relative coordinate n 2 is equal to / (^12) of the starting physical system, and 
is thus described by Eq. (fT5)) . 8 ' 9 For more than two electrons, in the case 
of a confined system (atom, molecule), the limit £ — > in Eq. p^|) describes a 
cluster of fermions, and Eq. (|15p becomes an approximation 8 ' 9 ' 17 for the internal 
degrees of freedom of the cluster. 

Here we focus on the high-density limit of the hamiltonians of Eq. ([5|) and 
we thus stick to the case N — 2. In general, the effective electron-electron 
interaction w c g(r 12) in Eq. (|15p is not known, and must be approximated. In 
the case of the He series, we have found 7-9 that a simple approximation based on 
the original idea of Overhauser 10 ' 11 gives very accurate results for 2 < Z < 10. 
In what follows we analyze the performance of the same approximation in the 
very Z — > 00 limit, and we extend our study to the k — > 00 limit of the Hooke's 
atom series. Of course, in the special case of the Hooke's series, the hamiltonian 
([5]) is exactly separable into center-of-mass and relative coordinates, so that 
the exact w e s(r '12) is directly available. However, the point here is to check 
whether the same approximate u> e ff(?"i2) that accurately describes the He series 
is capable to yield also good results for the Hooke's series, since this seems to 
be not the case for the currently available correlation energy functionals. 14 

The construction of an approximation for the e-e effective potential w e g 
starts with the decomposition 7-9 



w e ff(ri 2 ) = wfi(r 12 ) + < ff (ri 2 ), 



(17) 



where w£g 



V 2 v / /ks/V/ks is the potential that generates the Kohn-Sham 
/ks via Eq. (Tl5)) . and w^ s (ri 2 ) is a correlation potential that needs to be ap- 
proximated. In the usual DFT language, Eq. (fT7)) implies that we are treating 
exchange exaclty. 

In scaled units s, using standard perturbation theory we obtain the equation 
for the first-order contribution to / [see Eq. (jT2"j) ]. that separates into the Kohn- 
Sham and the correlation parts: 



^KS 



,(0) 



1>, 



(1) 



.(1) 
"KS ~ 

,(1) 



KS(1) 
°eff 



W 



c(l) 
eff 



^(0) 



(18) 
(19) 
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where ^ (0) = v 7 /^ /is = 2^ (0 Vks> = 2?/> (0 Vi 1 \ and 



„KS(o), o ^ _ 2(84 2 -8af 2 -38a; 2 -36ai 2 -9) 

(4 S 2 2 + 6s 12 + 3) 2 



= V " | , ,vr ~ -l (Heserie S )( 2 0) 

„2 



KS(0) (s 12 ) = (Hooke's series). (21) 



■a 



In Eq. ([IH]), /^g is exatly known for both series, so that we can also obtain 

KS(1) , . 

w cS by inversion. 

We thus concentrate on the correlation part, since we want to test approxi- 
mations for w° s . Defining u c (x) — xip^\x) and uq(x) = xip^(x), we have 



2 



jL _ w ks (°) + f (o) 

dx 2 " !i 



off 



wo- (22) 



Following the method of Refs. 15, 18, 19 we look for a solution of the kind 
u c {x) = uq{x) y(x). The function y{x) is then given by 

r x rl-r' C x ' 

V(*)= zSztJ ul{x")[w c c £\x")-e^]dx" + C 2 . (23) 



Iq u o( x ') Jo 

The constant C 2 is fixed by requiring the proper normalization, 

/>oo 

fi 1 \x)x 2 dx = => / y(x)ul(x)dx = 0. (24) 



The other integration constant has been fixed in Eq. ([23|) by setting equal to 
zero an unphysical term C\ f UQ 2 (x')dx' that would make u c (x) diverge for 
large x. 

3.2 Testing approximations: the Overhauser potential 

In Refs. 7-9 an approximation for w^g was built as an average "Overhauser- 
type" potential, 10 - 11 

«&(ri2) « f — + % ~ 4") * " ^12) - (25) 
Vri 2 2rf 2r s / 

where #(;c) is the Heaviside step function and r s is related to the average density, 
or, better to the dimension of the system. For two-electron atoms it was simply 
estimated as 7 

r s = (fn)- 1/3 , (26) 

n=l|drn(r) 2 . (27) 

The idea beyond this approximation is the following. The e-e correlation po- 
tential Wg ff (ri2) changes the Kohn-Sham / into the physical one, and must 



where 



G 




S 12 - Z r 12 



Figure 1: The correlation part of the intracule density, f c = f — / KS , divided by 
Z 2 , as a function of the scaled variable S12 = Zr 12 for the He isoelectronic series. 
The "exact" results are obtained from the accurate wavefunctions of Ref. 20. 
Approximate results at finite Z using the "Overhauser model" are taken from 
Ref. 7. The Z — 00 result corresponds to Eq. (|23|) with the potential of Eq. ([28]) . 



thus keep the information on the one-electron density (which is the same in the 
two systems) while turning on the e-e interaction I/V12. In Eqs. (f2l>l) - (f2"T)) this 
information is approximately kept via the average density n. 

In scaled units, the Overhauser potential to first order in a, to be used in 
Eq. ([19]). becomes 

ttfW « ( 7- + $k - w) 9 & - Sl2) ' (28) 

\si2 2 ^s s / 

where, if we adopt the prescription of Eqs. l[2"6 ]) -(|27 ]) . s s — 3 1 / 3 + 0(a) for the 
He series and s s — + O(a) for the Hooke's atom series. 

Equation (|23|) with the potential of Eq. ([28|) can be evaluated analytically 
as a function of S12 and s s for both series, although the final expressions are 
cumbersome and will not be reported here. The resulting fc for the He series 
is shown in Fig.[IJ together with the corresponding scaled quantity, Z~ 2 f c (s/Z), 
for some finite Z. [Since fc ■ (s) — lim^^oo Z f c (s), and/ c (s) = Z~ 3 f c (s/Z), the 
quantity to be compared with fc (s) is Z~ 2 f c (s/Z).] For finite Z we show both 
the "exact" result 20 and the approximate result obtained with the Overhauser- 
type potential of Eqs. (j^5]l - (f?7|) . We see that the Z dependence of the short- 
range part of f c is very well captured by this simple approximation. Figure [T] 
also suggests that the Z — > 00 limit of the short-range part of f c is well described 
by this approach. In Fig. [5] we show the result for /i 1 -* in the case of the Hooke's 
series from the Overhauser potential compared to the exact one, finding very 
accurate agreement. 

The KS and the correlation components of J*- 1 -* are shown in Fig. [3] for both 
series. We see that in the case of the He series the KS and the correlation 
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Hooke's series 




Figure 2: The correlation part of the first-order intracule, / c (si2) [see Eq. (|12p], 
for the Hooke's series. The exact values are compared with the results from the 
Overhauser-type approximation of Eq. ([28| . 




Figure 3: The decomposition of the first-order intracule intracule density 
/^(siz) [ see Eq. (fl2]l]: the Kohn-Sham part and the correlation part. 
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parts have roughly the same depth, while in the case of the Hooke's series the 
correlation part is much deeper than the KS one. This is due to the fact that the 
KS part gives the change in the e-e distance probability distribution only due 
to the first-order change in the one-electron density. In the case of the Hooke's 
series the first-order change in the density is much smaller, since the harmonic 
confining external potential is stronger than the Coulombic one. Indeed, the 
function x(s) of Eq. |8j in the case of the He series is about twice the one for 
the Hooke's atom series. 



4 Adiabatic connection and correlation energy 

The APD /i 1 ' ) (si2) gives the correlation contribution to second order to the 
expectation (V ee ) of the Coulomb electron-electron repulsion operator, V ee = 
l/ri2, 

(V ee ) = -1 [a(K e ) (1) +a 2 (Ke) (2) +0(a 3 )] , (29) 
or L J 

where (V ee )^ = (V ee )<gl + (V ee ) c 2) , and 



'KS 

\(2) 



{V et 



/ 47TS1 2 / C (1) (si2)dsi2. (30) 

Jo 



Our fc 1 * 1 from the Overhauser potential give (V ee ) c 2 ^ = —0.10256 Ha for the He 
sequence (to be compared with the exact 21 value, —0.09333 Ha), and (V ee )c = 
—0.10377 Ha for the Hooke's series (to be compared with the exact 16 value, 
—0.09941 Ha). The error is thus 9 mH for the He series and 4 mH for the 
Hooke's series. 

The correlation energy can then be otbained via the adiabatic connection 
formula of Eq. |T]), which for Ec reads 

where fc ^ is the first-order correlated part of / for the system with interaction 
a w x (si2). If we were able to calculate the exact fc for any u> A , the resulting 
E^ from Eq. (|31[) would be independent of the choice of w x . However, when 
we deal with approximations, we can obtain better results with some choices 
rather than others. 

As in Ref. 7, we build an Overhauser-type potential along the adiabatic 
connection as 

w^ x (s 12 ;s s ) = w x (s 12 ) -/ nw x (\s - s 12 |) ds, (32) 

J\s\<s e 

where, in scaled units, if we stick with the choice of Eqs. (|26p - ([27|) . n = (47r) _1 
for the He series and n = (47r 3 / 2 ) -1 for the Hooke's series. The idea behind 
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Eq. (J2U) is that the average density n (and thus the average s s ) is kept fixed to 
mimic the fact that the one-electron density does not change along the adiabatic 
connection while we turn on the e-e interaction. 

4.1 Linear adiabatic connection 

If we set w x (si2) = A/S12, Eq. (|32|) simply gives the Overhauser potential of 
Eq. (]28[) with a multiplying factor A in front. From Eq. (|23|) . we see that 

(2) (2) 

this corresponds to ^ = (V ee ) c /2, as in the exact case. I.e., the simple 
approximation of Eq. (I32p has the correct scaling behavior in the a — > limit. 

(2) 

Our result for Ec with the linear adiabatic connection thus gives an error of 
4.5 mH for the He series and 2 mH for the Hooke's series. 

4.2 The "erf" adiabatic connection 

A choice for w x that separates short- and long-range effects is the "erf" adiabatic 
connection, 7 ' 22-25 w x (si2) — erf(A Si2)/si2, for which Eq. ([3Tjl becomes 



The Overhauser-type potential corresponding to this interaction is reported in 
the appendix of Ref. 7. For the He isoelectronic series with 2 < Z < 10, the 
Overhauser-type approximation combined with the "erf" adiabatic connection 
gives 7 correlation energies with errors within 4 mH, better than the linear adi- 
abatic connection that gives errors within 10 mH. 

_____ r<2\ 

In the weakly-correlated limit, instead, we obtained, via Eq. (|3"3"|) , Ec = 

(2) 

—0.041 Ha for the He series and Ec — —0.046 Ha for the Hooke's series. The 
errors with respect to the exact values, 6 mH and 4 mH, respectively, are thus 
slightly worse than those obtained with the linear adiabatic connection. 

5 The LDA failure in the high-density limit: an 
analysis from the intracule density 

As a further element of comparison, we also computed the first-order /i (S12) 
within the local-density approximation (LDA), 



where g c (^i2', t s ) is the pair-correlation function of the uniform electron gas 26 of 
density n — (47rrg/3) _1 , and 




(33) 




(34) 




1-1/3 



k F (s) = [37r 2 n (0) (s)] 1/3 , 



(35) 
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Figure 4: The correlated part of the intracule density, / c {s\2 ) , of order a = 1/Z 
for the He series [see Eq. ([12]) ]: the present calculation is compared with the 
LDA approximation (panel a). Panel b shows the same quantities multiplied 
by 47TSi2: the integral under each curve gives the correlation part of the second 
order contribution to the expectation value (V ee ), which diverges in the case of 
LDA. 
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Figure 5: The correlated part of the intracule density, / c . (S12), of order a = 
fc -1 / 4 for the Hooke's atom series [see Eq. (fl2])]: the exact result is compared 
with the LDA approximation (panel a). Panel b shows the same quantities 
multiplied by Atts^: the integral under each curve gives the correlation part of 
the second order contribution to the expectation value (V ee ), which diverges in 
the case of LDA. 
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With these definitions, the density parameter r s of the uniform electron gas is 
locally proportional to a. We have numerically evaluated the right-hand-side 
of Eq. (|34p at smaller and smaller a (i.e., at larger and larger Z and k), for 
< s\2 < 5. As a decreases, the results tend to a well defined curve, shown in 
Figs. S] and together with the result from the Overhauser model (He series) 
and the exact result (Hooke's series). 

Since, as shown by Eq. (|34|) , the a — > limit corresponds to the r s — > limit 
of the uniform electron gas pair-correlation function g c , to better understand 
the LDA result for f c we now analyze more in detail the high-density behavior 
of g c . This analysis extends and completes the one done in Ref. 27. When 
r s — ► 0, the short-range part of g c scales as 

g c (x,r s — >0) = r s g'p (x) + 0(rf In r s ), x = r 12 /r s , (36) 

where the function g^ (x) does not depend explictily on r s and has been com- 
puted by Rassolov et al. 2S It is accurately recovered by the model g c of Ref. 26 
that we have used in the evaluation of Eq. (f34|) . The scaled variable x is locally 
proportional to the scaled variable s±2 [see Eq. (I3"5j) ]. Equation (|36|) thus shows 
that the short-range part (corresponding to values of the scaled variable x not 
too large) of g c in the r s — > limit has a scaling similar to the one of the He and 
Hooke's series in the a — > limit. This is also reflected by a good performance 
of LDA for S12 < 1, as shown by Figs. [4] and [5] 

However, the high-density electron gas is an extended system with important 
long-range correlations that are not present in finite systems like atoms and 
molecules. In fact, the scaling of Eq. (|3"6"|) is not valid when x ^> 1: it has been 
shown that the long-range part of g c scales as 26,29,30 

g c (x > l,r 3 ) -> r 2 s h(v), (37) 

where v is another scaled variable, v — ^/rlx 1 which is thus locally proportional 
to y/as\2- The function h(v) has the following asymptotic behaviors: 

« 1) oc ji* 2 , » 1) a w" 4 , (38) 

which are also correctly included in the model g c of Ref. 26. When r s — > 
(i.e., q -i 0), even for very large x the scaled variable v is small, so that the 
long-range (x ^> 1) behavior of g c is more and more dominated by the small v 
part of h(v), i.e., it behaves more and more like v~ 2 rather than like v~ 4 . It 
is this increasing dominance of the "short-range component of the long-range 
part" that causes the oc log(r s ) behavior in the correlation energy per electron 
of the high-density electron gas, and thus the divergence of the LDA correlation 
energy in the large-Z and large-A: limit of the He and Hooke's atom sequences 
(see, e.g., Ref. 31). In fact, when Z — > oo (or k — > oo), the high-density long- 
range behavior of g c affects the long-range part of / c LDA (si 2 ) in Eq. (M}. 

The small- 1> behavior oc v~ 2 of the function h(v) is related to the I/V12 di- 
vergence of the Coulomb potential at small r\2- For this reason, the oc log(r s ) 
high-density behavior of the correlation energy is still present in a uniform 
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electron gas with screened (or short-range only) Coulomb interaction (e.g., 32 
erfc(Ari2)/ri2), while is absent in an electron gas with long-range-only interac- 
tion (e.g., 33 ' 34 erf(Ari 2 )/ri 2 ). 

6 Conclusions 

We have computed the intracule density and the correlation energy for the 
high-density (weakly-correlated) limit of the He and Hooke's atom isoelectronic 
series via an approach 7-9, 17 based on an "average pair density functional the- 
ory" (APDFT), and inspired by the seminal work of Overhauser. 10-13 Unlike 
the currently available correlation energy functionals analyzed in Ref. 14, the 
APDFT approach gives accurate results for both series. In its present formu- 
lation, the APDFT approach works well for two-electron systems and for the 
uniform electron gas. Its extension to many-electron systems of nonuniform 
density is a big challenge, and we are presently exploring several different paths 
to achieve this ambitious goal. 9 

We have also analyzed the LDA failure in the same weakly-correlated limit of 
the He and Hooke's atom series, in terms of the long-range part of the intracule 
density. The results of Katriel et al. show that higher-order functionals such 
as PBE 35 and TPSS 36 can reasonably fix the LDA problems in the case of the He 
isoelectronic series, but are much less satisfactory for the Hooke's atom sequence, 
yielding a wrong scaling in the k —> oo limit (PBE) or a correct scaling with an 
error of about 40% on the asymptotic value of the correlation energy (TPSS). As 
stressed by Katriel et al., 14 these differences in performances for the two series 
raise serious doubts on the universality of currently available correlation energy 
functionals. The accuracy of the results obtained via the APDFT approach for 
both series suggests that the effort towards its generalization to many-electron 
systems of nonuniform density could be really worthwhile. 
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A fHs(rn) for the He isoelectronic series 

For the He isoelectronic series Eq. (fT3|) corresponds to 




Ei(-4x)] + e~ 2a: [164 + 27(3 + a;(9 + 4a;(3 + 2 a;)))[Ei(-2 x) 
7 - log(x)] + 3 x [-163 + 6 x (15 + x (7 + 10 x)) - 27 log(4/3)] 




1621og(2) + 811og(3)] , (39) 



14 



where 7 and the function Ei have been defined after Eq. ©. 
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